Probing Primordial Non-Gaussianity with Large-Scale Structure 



m 
o 
o 

(N 
o 

CD 

Q 



> 

oo 
(N 
(N 

o 
o 



Roman Scoccimarro"'^, Emiliano Sefusatti-'^ , Matias Zaldarriaga^ 
^ Center for Cosmology and Particle Physics 
Department of Physics, New York University 
New York, NY 10003 
^Department of Astronomy and Department of Physics 
Harvard University, 60 Garden Street 
Cambridge, MA 02138 

We consider primordial non-Gaussianity due to quadratic corrections in the gravitational potential 
parametrized by a non-linear coupling parameter /nl- We study constraints on /nl from measure- 
ments of the galaxy bispectrum in redshift surveys. Using estimates for idealized survey geometries 
of the 2dF and SDSS surveys and realistic ones from SDSS mock catalogs, we show that it is possible 
to probe |/nl| — 100, after marginalization over bias parameters. We apply our methods to the 
galaxy bispectrum measured from the PSCz survey, and obtain a 2(T-constraint |/nl| < 1800. We 
estimate that an all sky redshift survey up to z ~ 1 can probe |/nl| — 1- We also consider the use 
of cluster abundance to constrain /nl and find that in order to be sensitive to |/nl| — 100, cluster 
masses need to be determined with an accuracy of a few percent, assuming perfect knowledge of the 
mass function and cosmological parameters. 



I. INTRODUCTION 

Rapid progress in microwave background 
anisotropy experiments and large galaxy redshift 
surveys is providing high quahty data that can be 
used to test the nature of primordial fluctuations. 
The leading scenario for explaining the initial seeds 
for the formation of structure in the universe is 
inflation, a period of accelerated expansion in the 
early universe during which quantum fluctuations 
in a scalar field driving the expansion are stretched 
outside the Hubble radius and stay frozen until 
they cross back during matter domination and grow 
by gravitational instability. The predictions from 
this scenario have been worked out in great detail 
during the last twenty years, with most models of 
inflation predicting a scale-invariant spectrum of 
adiabatic Gaussian fluctuations. 

The Gaussian nature of primordial perturbations 
is a direct consequence of the slow-roll conditions 
on the inflation potential, required for the potential 
energy to dominate over the kinetic energy of the 
field and produce a sufficiently long period of accel- 
erated expansion. Under these circumstances, non- 
Gaussianities are very small, of the order of the tilt 
in the scalar spectrum 0|. The tilt is known to be 
rather small Q, |^ , in the language of Eq. Q below, 
/nl ^ 0.05. This bound can be relaxed to /nl 1 if 



higher dimensional operators are suppressed by the 
lowest possible scale consistent with slow- roll 0). 

There are several modifications to the basic 
physics of inflation that can lead to larger non- 
Gaussianities. In some set of these models there 
are additional light degrees of freedom during in- 
flation. For example, if the inflaton field has more 
than one component, it is possible to generate sig- 
nificant non-Gaussianity in the adiabatic component 
through the coupling to isocurvature components, 
see e.g. [EIEH- ^ second scalar field, usually called 
the "curvaton" could be light during inflation and 
come to dominate the energy density of the universe 
after the end of inflation before producing effectively 
a second reheating |i,|l[l3[llElll3- The density 
fluctuations we observe today could be due to fluc- 
tuations in the curvaton and could be non-Gaussian. 

In addition, recent work 0, 0, 0| suggested a 
new possibility in which the fluctuations are gen- 
erated during the reheating period when the in- 
flaton energy density is converted into standard 
model particles with a fluctuating decay rate. In 
all these models non-Gaussianities are primarily as 
given by Eq. Q and could naturally be of order 
/nl ~ 5 — 30 (151, Primordial non-Gaussianity 
at this level should be detectable through measure- 
ments of the bispectrum of the CMB [l^j llH and 
in galaxy surveys, as we shall show in this paper. 
Finally there are models in which the inflaton is not 



a slowly rolling scalar field but rather a fast moving 
ghost condensate [l^. In this case non-Gaussianities 
are much larger, close to the current upper limit but 
are not as simple as Eq. 

In this paper we consider departures of Gaussian- 
ity where the primordial gravitational potential at 
sub-horizon scales has the form |23| 

1>P™-'^+%(0'-(0')) (1) 

where (f) is a random Gaussian field, c is the speed 
of light, and for simplicity we assume /nl is a num- 
ber independent of scale. These are generically pre- 
dicted by all models in which the non-Gaussianities 
are generated outside the horizon |T5|. 

Measurements of the microwave background 
anisotropy bispectrum give la limits from COBE 
I/nlI < 1500 113, MAXIMA |/nl| < 950 |2^, and 
recently 2cr limits from WMAP of -58 < /nl < 
134 [is]. Upon completion, WMAP is expected to 
reach la sensitivity of order |/nl| ^ 20, whereas the 
Planck satellite would yield |/nl| ^ 5 |17l |. 

An alternative way of constraining primordial 
non-Gaussianity is by measuring the bispectrum of 
the galaxy distribution and looking for deviations 
from the predictions of gravitational instability from 
Gaussian initial conditions jl^, 13, (see e.g. [^^ 
for a review). In this case, there has been no con- 
straint yet on the particular model given by Eq. , 
although estimates have been made in |2fi| regard- 
ing the ability of large-scale structure to constrain 
/nl, concluding that galaxy surveys such as 2dF and 
SDSS will be able to probe only |/nl| ~ 10^ - 10^ 
and that galaxy surveys in general will not be com- 
petitive with CMB experiments in probing this type 
of non-Gaussianity. 

In this paper, we revisit the issue of how well large- 
scale structure can constraint non-Gaussianity of the 
type given by Eq. (Q, and reach a quite different 
conclusion. We show in particular that the SDSS 
galaxy bispectrum should be able to probe values of 
order |/nl| ^ 10^, and that a hypothetical all-sky 
survey with similar density up to z ~ 1 should be 
able to probe |/nl| ~ 1- We illustrate our results 
by applying these ideas to the measurement of the 
galaxy bispectrum in the PSCz survey (2^ . and find 
2(7 constraints |/nl| ^ 1800, comparable to the lim- 
its from CMB measurements before WMAP. 



The difference between our results and those 
in [2II can be traced to the assumption made in 
that constraints on primordial non-Gaussianity can 
be effectively "read off' from constraints on the non- 
linear bias parameter, which is independent of scale. 
This ignores the anomalous scale dependence of the 
bispectrum induced by primordial non-Gaussianity, 
which plays a crucial role in obtaining limits on 
it [2II23. 

Constraints on primordial non-Gaussianity of the 
type given by Eq. Q have been considered also by 
using gravitational lensing (29j . where it was found 

— 1/2 

that it is possible to achieve /nl ~ 150/gj^y , with 
/sky the fraction of sky covered, using lensing to- 
mography with 4 redshift bins up to ^ = 500. In 
addition, the use of the abundance of massive clus- 
ters to constrain /nl has been studied in [sollsH . In 
this work we also consider how well it is necessary 
to determine cluster masses to be able to use cluster 
abundance to probe /nl to the accuracy required by 
present upper limits. 

This paper is organized as follows. In sectionlTllwe 
discuss how the large-scale structure (LSS) bispec- 
trum is modified due to primordial non-Gaussianity 
given by Eq. Section Hill presents a signal-to- 
noise analysis for determining bias parameters and 
/nl from surveys, including a somewhat detailed cal- 
culation for the particular case of the SDSS survey, 
and application to the bispectrum of galaxies in the 
PSCz survey (Section UnEll. Finally, in Section HVI 
we consider the use of the abundance of clusters to 
constrain primordial non-Gaussianity. 

II. THE LSS BISPECTRUM 

Well after the universe becomes matter domi- 
nated, the fluctuations in the gravitational potential 
at time given by the scale factor a are related to the 
primordial fluctuations by 

'i>k(a) = ^^Tik) 1>r\ (2) 

where T(fc) is the transfer function, D^{a) is the 
growth factor linear perturbation theory, and $p''™ 
denotes the primordial gravitational potential at 
sub-horizon scales before matter-radiation transi- 
tion, and we have neglected anisotropic stresses. The 
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matter density is related to the potential by Pois- 
son's equation, which in Fourier space reads 



(5k(a) 



2 a A:2 



3 ^rriHc 



r$k(a), 



(3) 



where Jim is the present value of the dark matter 
density in terms of the critical density and _ffo = 
100 /iMpc~'^km/s is the present value of the Hub- 
ble constant. We assume cosmological parameters 
consistent with those determined by the WMAP ex- 
periment ,2], Vim = 0.27, f^b/i^ = 0.0224, h = 0.71, 
assuming a flat universe with a cosmological con- 
stant. The transfer function is computed using 
CMBFAST [s^, leading to a power spectrum nor- 
malization cTg = 0.82. Introducing 



M{k,a) = :j d^ 



5 n. 



(a) 



we can write 



(5k(a) 



M{k,a) ^f^, 



and in general 



(^ki5k....<5k«) = M{ki)M{k2)...M{kN) X 



V ki k2 "■ kjv ' 



(4) 



(5) 



(6) 



Henceforth we shall suppress the dependence on the 
scale factor, assuming a — 1, and drop the super- 
script denoting the primordial gravitational poten- 
tial, which is understood in all our expressions that 
follow. In our convention the power spectrum and 
the bispectrum are given by 

{SkM = SD{k,2) P{ki) (7) 
((5ki(5k2^k3> = 30 0^123) B{ki,k2,k3) (8) 

where k^.. j = k.^ + . . . -I- kj. We can write the lin- 
ear power spectrum of the density field as P^(fc) — 
M'^{k)P,j,{k) where P$(fc) is the primordial gravita- 
tional potential power spectrum, given by Eq. ^ 

P^{k) = p^{k) + 2^ J d'qP4q)P4\k~q\) 



(9) 



Equation |^ can be seen as a quadratic approx- 
imation to a more general power series expan- 
sion, for consistency we neglect higher-order correc- 
tions than those of leading order in the primordial 



non-Gaussianity parameter /nl- For example, the 
second-order correction in Eq. ^ changes the value 
of (Tg by less than 1% for /nl = 10^. From Eq. ||2J), 
it follows that the bispectrum in linear perturbation 
theory is given by 



^123 



M{ki)M{k2)M{k3) B^{ki,k2,k3) 
where -B123 = B{ki, k2,k3) and 
2/nl 



(10) 



B^{ki,k2,k3 



[P^{ki)P4k2) + eye: 



+o{fL)- 



(11) 



To get a sense of how significant is primordial non- 
Gaussianity for the density field as a function of /nl, 
we calculate the dimensionless skewness parameter, 
defined as 



(4)3/= • 

where the smoothed density field Sr is given by 



(12) 



Sr = J d^k W{kR) (5k, 



(13) 



where W{kR) denotes the Fourier transform of a 
spherical top-hat window of radius R in real space. 
From Eqs. H5I11I12|I we have 



S3{R) 



/NL 



{S'r 



2 \-3/2 



XM1M2A/12 W1W2W12 



Jd^ki d'k2P^{ki)P^{k2) 

(14) 



where = M{ki) and = W{hR) and ki2 = 
|ki -|- k2|. We can integrate this equation numeri- 
cally (30ii] , but it is also possible to derive an analytic 
expression that is exact at large scales and illustrates 
the basic dependence of S3{R) on cosmological pa- 
rameters. The non-trivial part in Eq. (|14|l is the 
integration over the angle between ki and k2, due 
to the dependence of M12 on the transfer function 
T{ki2)- We use the approximation. 



47r 



W12 fc?2 Tl2 



kfTiWi ( 1^2 
+ 1 ^ 2, 



(15) 



which is exact at large scales where T{k) is inde- 
pendent of fc, in view of the summation theorem of 



3 



0.05 




50 100 
R [Mpc/h] 



1000 



FIG. 1: The dimensionless skewness parameter S3 (J?) 
against smoothing scale R for /nl = —100 (solid line) 
and the approximation given by Eq. 1161 (dashed line). 



Bessel functions (see Appendix C in 26]). After sim- 
ple algebra, Eq. lfT3|l leads to 



/nl {Sr4>r) 



^ 1 d\n(SR(t)R} 
6 d\nR 



(16) 



This formula illustrates that the level of non- 
Gaussianity is proportional to /nl times the am- 
plitude of potential fluctuations smoothed on scale 
R, with a constant of proportionality that depends 
on the shape of the density-potential power spec- 
trum. Figure shows a plot of S3 (R) as a func- 
tion of scale R calculated by numerically integrating 
Eq. 114|) (solid line) and the analytic expression in 
Eq. H16I) . which matches the exact result at large 
scales. 

At the scales relevant for galaxy surveys, we must 
take into account second-order corrections in pertur- 
bation theory (PT), which read 



J d3gF2(q,k~q) 



(17) 



where (x = ki • ^2) 

F2(k„k2) = - + -(^ + ^)+-.^ (18) 

This gives for the power spectrum, 

P{k) = P^(fc) + 2 J Ai^2(k + q,-q)S^(k,q) 

= P^(fc)+P^(fc), (19) 

however, this correction is basically negligible at the 
scales we are interested, see Eq. (|^ below. For the 
bispectrum we have three contributions j23l | , 



(20) 



of which is the linearly evolved bispectrum due 
to primordial fluctuations, B"^ is the contribution 
due to the initial trispectrum, and B^ is the usual 
bispectrum generated by gravitational instability 
from Gaussian initial conditions. They are given by 



5^23 = 2F2(ki,k2)PiF2+ eye. 



(21) 



B 



123 



dV2(ki2 -q,q)T^(ki,k2,ki2 -q,q) 
eye. (22) 



where Pi = P(ki). The linearly evolved initial 
trispectrum is a quantity of second order in /nl 
and can therefore be ignored as we now show. We 
can estimate the corrections due to primordial non- 
Gaussianity for an equilateral triangle of side k (at 
large enough scales), 



B^ 

pB 
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e B"^ 

t 

-Se^ k^al B^ 



(23) 
(24) 
(25) 



where is the one-dimensional velocity disper- 
sion in units of the Hubble constant. 



q 



40 (/i"iMpc)2, (26) 



and the small parameter e is given by 
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^ ~ M{k)c^ ~ 
which leads to 



-0.07 



NL 



1 



h 

100 T(fc) 



0.01 /iMpc"^ 



ek^al ~ 3 X 10" 



/nl 
100 



T{ky 



(27) 



(28) 



Equation (|25|l says that the power spectrum cor- 
rection due to primordial non-Gaussianity, Eq. H19() , 
is suppressed by e k^al at large scales, and becomes 
of order a few percent as fc > 0.1 /iMpc~^, inducing 
a scale dependence on the bispectrum, see Fig. |3 
below. At these scales other effects due to non- 
linearities (specially redshift distortions) become im- 
portant, we found that in redshift space it is very 
difficult to see the effects of nonzero /nl at these 
scales. Equation (|24|l shows the trispectrum cor- 
rection given by Eq. H22|l is negligible. Therefore 
we only probe primordial non-Gaussianity through 
Eq. 

It is convenient to introduce the reduced bispec- 
trum llll, defined as 



Q 



23 



123 



P1P2 + P1P3 + P2P3 



(29) 



which is independent of time for Gaussian initial 
conditions; moreover, in tree-level perturbation the- 
ory it reduces to a scale independent value for equi- 
lateral configurations |33J 



(fc) 



B^{k,k,k) 

T[p(fcjF 



(30) 



For general triangles, Q retains approximately this 
simple behavior, it is independent of power spectrum 
normalization, and only very weakly dependent on 
f2m through the factor ~ Q,m^^^ the only rele- 
vant dependence of Q is on the local spectral index 
neff(fc) = dlnP/dlnfc and triangle shape through 
Eq. 

An important consequence of primordial non- 
Gaussianity from Eq. ^ is that it violates the 
scaling induced by gravity, since Q^qik) ~ e ~ 
l/[fc2r(fc)], see Eq. The top panel in Fig. H il- 

lustrates the deviations from Eq. H30|) when |/nl| = 



10^. Note that the scale dependence seen here is op- 
posite to that in the skewness (compare to Fig. 
this is simply due to the difference in normaliza- 
tions between Eq. H12|l and (|29|l . The bottom panel 
in Fig. 121 shows the corrections due to primordial 
non-Gaussianity as a function of triangle shape for 
fci = 0.02 /i Mpc""'^ and fc2 = 2fci. We now explore 
how well these deviations can be probed with galaxy 
surveys. 



III. BISPECTRUM SIGNAL TO NOISE 

A. Order of Magnitude Estimate 

In this section we will calculate the minimum /nl 
that can be measured by a galaxy redshift survey as 
a function of the survey parameters. We start by 
making a simple estimate to put the results of the 
next sections in context. A survey with volume V 
contains 



An 



k-^ 



V 



(27r)3 



4.5 X 10^ 



V 



(l/i-iGpc)3 V0.3/iMpc 



(31) 



Fourier modes, where fcmax is the largest wavenum- 
ber that can be used in the analysis. Let us use 
the skewness to estimate /nl- With this number of 
modes we expect to be able to measure the skewness 
roughly to the level, 

/ 1 c; 

(32) 

where we have used the Gaussian variance. In the 
previous section we showed that the skewness is of 
order S3 2 x 10~'*/nl so we expect to be able to 
detect /nl ~ 30 ff F ~ (1 h-^Gpcf. 

The above estimate indicates that we expect LSS 
surveys to be competitive with CMB experiments in 
constraining primordial non-Gaussianity. However, 
there are several simplifications in this estimate: 1) 
The density field fiuctuations are not Poisson dis- 
tributed, thus the skewness is not the best estimate 
of /nl; 2) The galaxy density field is a biased tracer 
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k [h/Mpc] 



of the underlying mass and biasing modifies the bis- 
pectrum, therefore one must determine simultane- 
ously /nl and bias parameters; 3) The survey geom- 
etry can significantly complicate the determination 
of the bispectrum. We will tackle these problems 
in the rest of this section to obtain a more robust 
estimate of the capabilities of redshift surveys. 



B. Ideal Geometry 

We first consider the case of ideal survey geometry, 
assuming that bispectra for different triangle shapes 
are uncorrelated, i.e. the bispectrum covariance ma- 
trix is diagonal and given by Gaussian statistics. In 
order to see how well one can probe non-Gaussianity 
one has to include also the possibility that galax- 
ies are biased tracers of the density field. At large 
scales, it is reasonable to assume that biasing is lo- 
cal, then Isi 



1 . r 

1.6 


1 1 


1 


1 1 1 1 
k|=0.02 h/Mpe 


1.5 


\ \ 




k3=k3{k,,k,,e) ^/ /I 

' / 


1.4 


^ \ \ 




1 / / 
/ / / - 


1.3 






/ / / 
// / 
/ / / 


1.2 


% 




/ / / 
II 
II 
II 


1.1 






II 
II 
II 


1 


\\ 
\ 

\ 


\ 

\\ 


' 1 
' 1 


0.9 
0.8 
0.7 




\ \ 

\ 

\\ ^ 
\\\ 


' 1 1 

'/' 
'11 

/ /' 

■ / / 


0.6 

( 




1 


""i' 1 1 1 


) 0.1 0.2 0.3 


0.4 0.5 


0.6 0.7 0.8 0.9 



e/7i 

FIG. 2: Primordial Non-Gaussianity corrections to the 
reduced bispectrum for equilateral configurations as a 
function of the wavenumber k (top panel) and for k\ — 
0.02/iMpc~^ and = 2fci as a function of the angle 6 
between ki and k2 (bottom). 



Sg=h6 + ^S^ + ... (33) 

where bi and 62 are constants. The bispectrum in 
the galaxy distribution, including primordial non- 
Gaussianity, will be given by 

Bg{ki,k2, fcs) = bl Bg^ + blb2 (F1P2 + eye.) + 
+bl B^23 (34) 

In terms of the galaxy reduced bispectrum Qg, we 
have 

Qg~ +^2+ f,^ ■ (35) 

Note that each term in this expression has a different 
behavior, Qf23 depends very weakly on scale through 
the local spectral index (which can be measured) and 
depends strongly on triangle configuration, the sec- 
ond term due to non-linear bias is a constant, and 
the last term due to primordial non-Gaussianity de- 
pends rather strongly on scale (see Fig. [21). Therefore 
it is possible to simultaneously obtain constraints on 
61, &2 and /nl- 

For the reasons discussed above, we work with 
the reduced bispectrum Q123 |6l| . which has identi- 
cal signal to noise to B123 in the limit of Gaussian 
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fluctuations. Indeed, for Q123 = -B123/S123, where 
E123 is the denominator in Eq. {AQ'^)/Q'^ = 

(AB2)/B2 + (AE2)/S2, and {AB^)/B'^ ~ [3A(fc)]-i 
[see Eq. ^ below] and (AI]^}/^^ ~ where 
Nk is the number of k-modes contributing to the 
estimate of the power spectrum P{k) and A(A:) = 
ATTk^P{k). In other words, the signal to noise of Q 
is dominated by that of B and the power spectrum 
can be considered perfectly determined for our pur- 
poses. 

The bispectrum signal to noise for a given triangle 
can be written as, 

f —\ = Q123 ^ B123 , , 

^7Vyi23" AQ123 " AB123' ^ ' 

where the last equality follows from the discussion 
above. The bispectrum variance in the Gaussian 
limit can be computed in similar fashion to the stan- 
dard power spectrum case 36]. For a bispectrum 
estimator [87j 

B123 = -r^ [dPqi U^q2 [d^qs Sd {(1123) 

V123 Jki Jk2 JH 

X (5qi<5q2^q3 (37) 

where the integration is over the bin defined by qi G 
[ki - 5k/ 2, ki + Sk/2), Vf = {2tt)^/V is the volume 
of the fundamental cell, and 
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FIG. 3: Minimum detectable value of /nl as a function 
of fcniax (left) and fcmin (right), after marginalization over 
bias parameters. The solid line corresponds to an ideal 
survey with V = 1 (Gpc/h)^ and no shot noise. The 
other lines correspond to V = l{h~^ Gpc)^ with n = 
10-*(/iMpc-^)^ (dotted line), V = 0.3(/i-^ Gpc)=' with 
n = 3 X 10"^(/iMpc"^)^ (long dashed line) and V = 
0.1(/i"^Gpc)^ with n = 10"^( /iMpc-^)^ (short dashed 
line) . 



Vi23 = d^qi d^q2 d^qs Sd{cIi23) 

J ki ^2 ks 

~ Stt^ fcifc2fc3 Sk^, (38) 
the variance is 

AB^ ^ fc3 ^ Ptot{ki)Ptot{k2)Ptotik3) (39) 

1^123 

where S123 = 6, 2, 1 for equilateral, isosceles and gen- 
eral triangles, respectively, and 

P,„t(fc)^P(fc) + -i^ i (40) 

where the number density h accounts for the shot 
noise. 



For simplicity, in this section we shall assume that 
the estimates of the bispectrum are Gaussian dis- 
tributed (in practice one can check this assumption 
for a given survey geometry, see ) , we shall go be- 
yond this in the next section when we obtain bounds 
on /nl from the PSCz survey. In the Gaussian ap- 
proximation, the likelihood for the bispectrum esti- 
mates C obeys |3^ 

m r J- I \ ^ (Qobs Qmod) f at\ 

- 2 In £ = const + 2^ , (41) 

where T denotes sum over triangles [defined pre- 
cisely below, Eq. l|IS|l ]. Qobs is the observed bis- 
pectrum and Qmod and AQ^^^ are computed from 
Eqs. H35|) and H39|l in terms of the model parameters 
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bi, 62, /nl- If observations are consistent with the 
fiducial model with 61 = 1, &2 = and /nl = 0, 
Qobs — and we have 



— 21n£ = const + aiUj Fij, (42) 

where ai = (1 — &i), a2 = 62/^1, ^3 = /nl and the 
Fisher matrix is given by 



femax 

E 



E 



fe2 

E 



-"123^123 ^^2^ 
^-3^23 



with k* 



- max(A:niin, \k\ — ki\), and we have as- 
sumed that the variance A_B^23 is computed only in- 
cluding linear bias. Here -8^23 = ^1231 -^123 — ^123, 
and B 



(3) 
123 



^123 //nl, 



which are respectively the 
bispectra induced by gravity, non-linear bias and pri- 
mordial non-Gaussianity. Equation H42II is now the 
standard Gaussian likelihood for the parameters ai 
and their error bars (marginalized over all other a-,) 
are simply given by 



(44) 



Figure shows the minimal detectable value of 
/nl, given by Itr error bars from Eq. H44|) . for dif- 
ferent survey volumes and galaxy number densities. 
The left panel shows how the /nl limits improve as 
we include more triangles towards smaller scales by 
increasing fcmax, whereas the right panel shows the 
opposite regime, where fcmax is held fixed and one 
probes larger scales (from right to left) by decreas- 
ing fcniin all the way up to the fundamental mode 
of the survey fcmin = ^/ • Figures 13 and 01 assume 
bispectra whose sides are binned with 6h — kf. 

The different lines in Fig. |2| have been chosen to 
roughly represent the 2dF survey (short-dashed), 
the main sample (long dashed) and LRG sample 
(dotted) of the SDSS survey, and a hypothetical 
survey with the same volume as the LRG sample 
but with high enough density to make shot noise 
negligible at kj^a^- In the absence of shot noise 
(and keeping our ideal survey geometry constant) 
the minimum detectable value of /nl scales simply 



as y^^/^. The scaling with /cmax is basically given 
by the naive expectation that the constrains on /nl 
should be inversely proportional to the square root 
of the number of modes available, Nk cx k'^^^V. 
We thus see from Fig. O that an all-sky survey with 
n ~ 3 X 10~^ (/iMpc~^)^ up to redshift z ^ 1 can 
probe values of /nl of order unity. A redshift sur- 
vey of such a volume may be realistic in the not too 
distant future li^. 

Figure 0] shows the corresponding results for the 
bias parameters bi and 62, showing that 2dF and 
SDSS surveys should be able to determine the bias 
parameters to within 1-2 percent accuracy. This is 
in rough agreement with previous estimates |39l l4lj | , 
we will provide a more detailed assessment for the 
SDSS case in Section ITTTn 

It is worth comparing the right panels of Figs. 13 
and 21 to see from what scales is most of the signal 
coming from. For the case of biasing parameters, 
the dependence on fcmin at large scales is rather weak 
compared to that for /nl, saying that large-scale tri- 
angles contribute significantly more information to- 
ward constraining primordial non-Gaussianity than 
biasing parameters. This is again a consequence of 
the scale-dependence of primordial non-Gaussianity, 
in fact the signal to noise for an equilateral triangle 
of size k is given by (A = Ank^P) 



B'^ 4 ; 



(45) 
(46) 



for linear and quadratic bias, respectively, whereas 
for primordial non-Gaussianity, 



B^ 
AB 



ASiT^k^Ps 



/nl = 1.4 X 10 ''/nl = const, 

(47) 

where we have used a scale invariant primordial spec- 
trum normalized to produce the correct level of CMB 
anisotropics, {Airk^P^/c^y^^ ~ 2.3 x 10~^ Equa- 
tion l|47|) says that the signal-to-noise per triangle 
is constant 6^. In other words, for primordial non- 
Gaussianity given by Eq. ^ the total signal to noise 
is only decreased at large scales by the decline in 
the number of triangles, whereas for biasing param- 
eters there is an additional suppression due to the 
decrease in the rms fluctuation amplitude at large 
scales laSEl. 
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FIG. 4: Same as Fig. |3 but for the linear and 
quadratic bias parameters (62). 



TABLE I: SDSS moclc catalogs (for each geometry) and 
bispectrum measurements used in the analysis. Cosmo- 
logical parameters are as in Sect.|n] fcmin = 0.02 h Mpc~^ 
and fcmax ~ 0.3 Mpc~^ . 



10-Vnl 


mocks 


10-3Po 


-^triangles 


Nt 


0,1,4 


6080 


2,5,10 


7.5 X 10^° 


1015 



This explains why our constraints on /nl are 
about two-orders of magnitude better than those ob- 
tained in '25l , where the bound on /nl is derived by 
translating the constraint on 62 to an effective value 



of /nl at fc ~ 0.6 h Mpc 



C. Including Survey Geometry: SDSS Forecast 

Let us improve the above treatment, considering 
a realistic survey geometry with the induced covari- 
ance matrix between different triangles. We also in- 
clude redshift distortions, as calculated by second- 
order Lagrangian Perturbation Theory (2LPT) from 
non-Gaussian initial conditions given by Eq. l(T]l. 
see for a comparison of 2LPT against N-body 
simulations for the redshift-space bispectrum. For 
biasing, we assume Eq. H35() still holds in redshift 
space, which is a reasonable approximation near our 
fiducial unbiased model. A treatment of bias and 
primordial non-Gaussianity in the presence of red- 
shift distortions is beyond the scope of this paper. 

We consider two survey geometries that approx- 
imate the north part of the SDSS survey, a 7, 300 
square degree region IgII and a second one with 
10,400 square degrees [63 • We don't include the 
South part of the survey in our analysis, which has 
a smaller volume and a nearly two-dimensional ge- 
ometry that complicates the simplified bispectrum 
analysis we will do below. For the radial selection 
function we use that corresponding to the "NYU LSS 
Samples" 10-12 and we assume that the angu- 
lar selection function is unity everywhere inside the 
survey region, which is a very good approximation. 

Using a 2LPT code [H with about 42 x 10^ 
particles in a rectangular box of sides Li = 
660,990, and 1320 /i"""^ Mpc, we have created about 
6 X 10^ realizations of each geometry, for Gaussian 
initial conditions and models with primordial non- 
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Gaussianity with /nl — 100 and 400 (see Tabled). 
In all cases, cosmological parameters are as given in 
Sect .mi and bi = 1, 62 = 0. For each of these realiza- 
tions, we have measured the redshift-space bispec- 
trum for triangles of all shapes with sides between 
fcmin = 0.02 /iMpc^^ and fcmax = 0.3/iMpc~^, giv- 
ing a total of -/Vtriangics = 7.5 X 10^° triangles. These 
are binned into A^t = 1015 triangles with a bin 
size of 5k = 0.015 /iMpc~^. The generation of each 
mock catalog takes about 15 minutes, and has about 
4 X 10^ galaxies for the smaller area and 5.7 x 10^ 
galaxies in the larger area case. The redshift-space 
density field in each mock catalog is then weighed 
using the FKP procedure 0|, see e.g. [H, l39| for 
a discussion in the bispectrum case. We have tried 
different weights Pq (see Table to minimize the 
error bars, the results we present correspond to 
Pq = 5000 {h-^Mpcf. The bispectrum in each 
realization is then measured for all (~ 7.5 x 10^" 
binned into 1015) triangles, taking about 2 minutes 
per realization |6q . 

In order to generalize the discussion given above to 
the case of arbitrary survey geometry, we introduce 
the bispectrum eigenmodes g„ f38| . 

EQm ~ Qm , 
7m«— TT^ , (48) 

m=l 

where = ((3,„), (AQ„)2 = ((g„j - Q„,Y). By 
definition they diagonalizc the bispectrum covari- 
ance matrix, 

{q,i qrn) = Srun, (49) 

and have signal to noise. 




The eigenmodes are easy to interpret when or- 
dered in terms of their signal to noise [s^l- The 
best eigenmode (highest signal to noise), say n = 1, 
corresponds to all weights 7^1 > 0; that is, it 
represents the overall amplitude of the bispectrum 
averaged over all triangles. The next eigenmode, 
n = 2, has 7,^2 > for nearly coUinear triangles 
and 7m2 < for nearly equilateral triangles, thus it 



represents the dependence of the bispectrum on tri- 
angle shape (see Fig. (SJ. Higher-order eigenmodes 
contain further information such as dependence of Q 
with scale, important to constrain primordial non- 
Gaussianity [23, 123 1 see Fig. below for illustration 
of this point. 

If the bispectrum likelihood is Gaussian, we can 
write down the likelihood as a function of the pa- 
rameters aj as, 

Ci{a,}) PmiaM (51) 

1=1 

where the Pi (x) are all equal and Gaussian with unit 
variance. We have checked from our Monte Carlo 
pool that the distribution of Q is indeed Gaussian 
even at the largest scales we consider. In practice, 
Gaussianity of Pi is not guaranteed at large scales 
due to the deviations from the central limit theorem 
by lack of enough uncorrelated triangles 28, 38] . If 
not Gaussian, diagonalization of covariance matrix 
does not guarantee independence of the eigenmodes, 
thus C does not necessarily factorize as in Eq. if^ . 
but this is a good approximation for small deviations 
from Gaussianity when the non-Gaussian Pi{x) are 
determined from mock catalogs 13^ . 

We calculate the bispectrum x covariance 
matrix from our A'mocks realizations of the survey 
(see Table and from that obtain 7m„ and A„, 
which gives the ingredients to implement Eq. H51|l . 
The results from such likelihood analysis are shown 
in Fig. [S] Contours denote joint 68% probabilities, 
two parameters at a time marginalized over the third 
parameter. The inner contour corresponds to the 
survey geometry with the larger area. The lower 
right panel shows the resulting la error bars after 
marginalization, smaller uncertainties correspond to 
the larger volume survey geometry. We have scaled 
our /nl = 100 bispectrum measurements to contin- 
uous values of /nl, identical results are obtained by 
scaling the /nl — 400 mock catalogs. 

It is difficult to compare these results to those of 
the previous section, since they correspond to very 
different survey geometries. However, comparing 
Fig.Elto the long-dashed lines in Figs. l3l4l shows that 
our more realistic estimates give error bars larger by 
a factor of 4-5. There are reasons to expect our 
"realistic" treatment to be actually an upper bound 
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FIG. 5: Joint 68% confidence intervals for two param- 
eters at a time, witli tlie tliird parameter marginalized 
over. The inner contour corresponds to the larger survey 
volume case. The lower right panel shows the resulting 
la error bars after marginalization, with smaller errors 
corresponding to the larger volume survey geometry. 



FIG. 6: 95% confidence limits on /nl from the PSCz 
galaxy bispectrum |27| after marginalization over bias 
parameters, as a function of the number of eigenmodes 
included in the likelihood analysis. 



to the achievable error bars with a more sophisti- 
cated analysis, for the following reasons. First, we 
only considered the north part of the survey; second 
we use FKP weighting, which is not optimal at large 
scales and thus could potentially reduce our sensitiv- 
ity, particularly to primordial non-Gaussianity; and 
finally, we have only used closed triangles in Fourier 
space, due to the lack of translation invariance there 
is also signal in open configurations. 

It is interesting to compare the results of Fig. O 
between the two geometries. The larger volume sur- 
vey leads to an improvement in marginalized error 
bars of 20% for h, 35% for &2 and 45% for /nl- 
This is more than what one expects for uncorrelated 
contributions to the constraining power of the sur- 
vey due to the increased volume (~ 20%), and is 
a manifestation in the improvement of the bispec- 
trum covariance matrix due to the narrower survey 
window function in Fourier space. 



D. /nl from the PSCz galaxy bispectrum 

We now consider constraints on /nl from the 
gal axy bispectrum measured in the PSCz survey 
in jig. The PSCz survey Q is based upon the 
IRAS Point Source Catalog, the bispectrum mea- 
surements we use are based on 13180 galaxies in the 
range 20 h^^ Mpc < i? < 500 /i"^ Mpc, with galactic 
latitude |6| > 10°, and IRAS 60 micron flux /go > 
0.6Jy, see for more details. We use triangles 
from fcmin = 0.05/iMpc~^ to fc^ax = 0.3/iMpc~"^, 
binned into A^t = 183 triangles. We use the eigen- 
modes and their probability distributions Pi (x) com- 
puted from ^ 10^ 2LPT mock catalogs in [38j . 

Figure El shows the 95% confidence limits on /nl 
after marginalization over the bias parameters, as 
a function of the number of eigenmodes (ranked by 
signal to noise) included in the likelihood analysis. 
Note how the constraint on /nl is set by the n > 2 
eigenmodes, which are sensitive to scale dependence 
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of the bispectrum. The bound on /nl converges 
after the best ~ 100 eigenmodes are mcluded, since 
the remaining half of the eigenmodes does not add 
any significant signal to noise. The 95% confidence 
limits corresponding to 123 eigenmodes are 



2000 < /nl < 1600 (95%). 



(52) 



This is comparable to the constraint from CMB fluc- 
tuations before WMAP although LSS is 
sensitive to somewhat smaller scales than the CMB. 
Our analysis in the previous section suggests that 
these limits should be improved by about an order 
of magnitude by the 2dF and SDSS surveys. 



IV. CLUSTER ABUNDANCE AND 
PRIMORDIAL NON- GAUSSIANITY 

The abundance of clusters probes the probability 
distribution function (PDF) of the matter fluctua- 
tions, and it is thus a natural candidate to constrain 
primordial non-Gaussianity |M ill E El El El ■ 
For large masses, the abundance of clusters depends 
on the right tail of the PDF which decays exponen- 
tially for Gaussian initial conditions. However, be- 
fore this can be used to place a constrain on /nl, 
it is necessary to have under control a number of 
systematic effects. 

Even with recent progress in the determination of 
cosmological parameters, uncertainties on 17^ and in 
particular erg can alter the Gaussian abundance pre- 
diction enough to make difficult probing the small 
levels of primordial non-Gaussianity corresponding 
to e.g. /nl — 100. In addition, calculation of the 
mass function from a given PDF is not straightfor- 
ward, even for Gaussian initial conditions measure- 
ments in numerical simulations suffer from system- 
atic uncertainties of order 10 — 30% depending on 
the definition of halo mass ^49,] . For primordial non- 
Gaussianity of the type given by Eq. (Q) , the mass 
function has been estimated analytically in |30l | , but 
there has been so far no complementary study using 
numerical simulations to give a sense of the uncer- 
tainties involved. 

At a more fundamental level, it is not even clear 
that one can probe /nl using rare events such as 
clusters, given that the tail of the PDF in general is 
not determined by the skewness alone, but rather 



all the higher-order cumulants of the distribution 
through its generating function p^ . In other words, 
the mass function for massive clusters will in gen- 
eral depend on higher-order than quadratic terms 
in Eq. On the other hand, the first step is to 
see whether deviations from Gaussianity can be de- 
tected at all for models such as Eq. (Q, keeping in 
mind that predictions for the mass function in a fully 
specified model $(</)) may be rather different than 
truncating Eq. to second order in 0. 

An important source of uncertainty in comparing 
theoretical predictions with observations of cluster 
abundance is the determination of cluster masses. 
Here we will ignore all the sources of uncertainty 
mentioned above, and estimate to what accuracy 
it is necessary to determine cluster masses in or- 
der to distinguish a Gaussian from a model with 
/nl = —50,-100, which being positively skewed 
gives a larger abundance than the Gaussian case. 
For simplicity we assume the Press-Schechter (PS) 
formula for the mass function Ho, E3l , 



dn{M,z) _ 2p 
dM ~ M2 



dlna-M(^) 



dlnM 



iyP{i^) (53) 



where the mass M = AnpR^/S is related to the 
smoothing scale R and mean density p by a top 
hat filter and v = dc/crM{z), with Sc — 1.686 and 
aM{z) the variance of the density field at redshift 
z smoothed at scale R. A more accurate analytic 
estimate of the mass function in the Gaussian case 
is given by the Sheth-Tormen (ST) m ass function 
based on ellipsoidal collapse |52,|53,|5J|, but it is un- 
clear how to generalize it for the non-Gaussian case, 
therefore we will use Eq. (|53|l instead by changing 
the PDF P{i'). This has been found to be a rea- 
sonable approximation when compared to numeri- 
cal simulations, at least for large in non-Gaussian 
models with initial conditions fs^. We will only 
deal with ratios of non-Gaussian to Gaussian abun- 
dances, therefore our estimates should be more ac- 
curate than differences in the absolute calibration of 
the mass function. When estimating the total num- 
ber of clusters expected in the Gaussian case, we will 
quote both PS and ST values. 

The total number of clusters with mass larger than 
M that can be observed in an all-sky survey between 



12 



redshifts z\ and zi is given by, 



where 



with 



dV _ cDl{z) 



dzdfl H„E{z) 



dz' 

Ho Jo W) 



(55) 



(56) 



and E{z) — ^ilA/(l + z)^ + Qa- We assume a flat 
Universe with the cosmological parameters given in 
Sect.ini 

Figure [7| shows the excess in the abundance of 
clusters due to non-Gaussianity, given by 



. _ iV,r (> M) 

^NGihL) = M) 



(57) 



for /nl = —100 (dashed) and /nl = —50 (dotted), 
compared to the excess that may come from an un- 
derestimate in the mass determination, given by 



Ag{AM/M) 



^ Ng{> Af - AM) 



(58) 



plotted as a function of /S.M/M. To use Eq. we 
have measured the non-Gaussian PDF for the den- 
sity field from realizations of Eq. iQJ on a 512'^ grid in 
a box of 1 Gpc a side, smoothing the density field 
on a scale of 10 Mpc (other smoothing lengths in 
the relevant range do not change the results). Our 
results in Fig.[3for the excess abundance are in rea- 
sonable agreement with Fig. 8 in [8(tI |. where an an- 
alytic approximation was developed to calculate the 
mass function instead of using Eq. (|53|) . 

We see from Fig. [7| that in order to probe pri- 
mordial non-Gaussianity at levels |/nl| ^ 100, un- 
certainties in the mass limit M should be well be- 
low 10%, assuming perfect knowledge of cosmologi- 
cal parameters, and negligible cosmic variance. At 
present, direct mass determinations through weak 
lensing seems to suffer an absolute uncertainty not 
below AM ~ lO^'*/i"^M0 and larger at high red- 
shift with a contribution of distant large scale 
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FIG. 7: The horizontal lines show the excess in the num- 
ber of clusters in an all-sky survey with redshift lim- 
its as shown for masses above lO" and lO^'^ /I'^M© 
expected for non- Gaussian primordial fluctuations with 
/nl = -100 (dashed) and /nl = -50 (dotted). The 
solid curves show the excess in the number of clusters in 
the Gaussian case due to an underestimate of the mass 
hmit M by an amount AM as a function A.M/M. The 
number of clusters in the Gaussian case No for such a 
survey is given assuming the PS and ST mass functions. 



structure of the order of the 6% alone |52J , while the 
estimates provided by X-ray temperature measure- 
ments have statistical errors of 10 — 30%. 

It is interesting to note here that future cluster 
surveys using the Sunyaev-Zel'dovich effect 58] re- 
quire similar accuracy in the determination of the 
mass threshold for detection. The goal for these ob- 
servations to be able to probe cosmological param- 
eters is to achieve an accuracy of 5% on the mass 
hmit and 10% on the mass function [5^|63- These 
turn out to be minimal requirements for the use of 
clusters to probe primordial non-Gaussianity at lev- 
els comparable to that within reach of CMB and LSS 
methods. 
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V. CONCLUSIONS 

We studied constraints on primordial non- 
Gaussianity, of the form given by Eq. from mea- 
surements of the bispectrum in galaxy redshift sur- 
veys. We find that taking into account the scale- 
dependence of the bispectrum induced by primor- 
dial non-Gaussianity is essential to obtain reliable 
constraints on /nl- As a preliminary application, 
we derived the first constraints on /nl from LSS 
using the galaxy bispectrum measured in the PSCz 
survey, obtaining -2000 < /nl < 1600 at the 95% 
confidence level. 

We estimate that the SDSS survey should achieve 
68%-level constraints for |/nl| of at least 150-200. 
The uncertainty in this number is due to the simplifi- 
cations used in our analysis. Although we have taken 
into account the geometry in detail by using mock 
catalogs with the same selection function, there are 
many features of our analysis that need to be im- 
proved for a more accurate assessment. We only 
consider the north part of the survey, since the south 
region has a geometry that is nearly two-dimensional 
and that would invalidate our analysis that assumes 
the window function of the survey is sufhciently nar- 
row in all directions. To make our analysis more 
tractable, we have used a weighting scheme that is 



not necessarily optimal at large scales, for primordial 
non-Gaussianity this can make a difference since sig- 
nificant information is coming from large scales. For 
these reasons, our estimate is likely to be an upper 
limit to the value of /nl to be probed by SDSS. 

In addition, we showed using simple signal to 
noise estimates that an all-sky survey with n ~ 
3 X lO^'^ ( /iMpc~^)^ up to redshift z ~ 1 can probe 
values of /nl of order unity. A redshift survey of 
such a volume may be realistic in the not too dis- 
tant future [ioj . 

We have also studied the use of cluster abundance 
to constrain primordial non-Gaussianity. At present, 
uncertainties appear too large to be able to compete 
with CMB or LSS, but this can change if cluster 
masses can be determined with an accuracy of a few 
percent, and other systematic errors in our theoret- 
ical understanding improve. 
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